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Absolutely Maximally Entangled (AME) states are those multipartite quantum states that carry 
absolute maximum entanglement in all possible partitions. AME states are known to play a relevant 
role in multipartite teleportation, in quantum secret sharing and they provide the basis novel tensor 
networks related to holography. We present alternative constructions of AME states and show their 
link with combinatorial designs. We also analyze a key property of AME, namely their relation to 
tensors that can be understood as unitary transformations in every of its bi-partitions. We call this 
property multi-unitarity. 


I. INTRODUCTION 


A complete characterization, classification and 
quantification of entanglement for quantum states re¬ 
mains an unfinished long-term goal in Quantum Infor¬ 
mation theory. Nevertheless, a large number of rele¬ 
vant results related to entanglement are well-known by 
now. It is, for instance, known that generic condensed- 
matter physical systems characterized by homoge¬ 
neous nearest neighbor interactions carry an amount 
of entropy bounded by an area law. This fact opens 
the possibility of applying tensor network techniques 
to describe the ground state of relevant physical sys¬ 
tems, including such phenomena as quantum phase 
transitions. 

However, there are many situations where entangle¬ 
ment entropy is maximal or nearly maximal, that is, 
it scales as the volume of the quantum system. This 
is the case of generic time evolution of local Hamil¬ 
tonians, of evolution of a quantum computer address¬ 
ing a quantum Merlin-Arthur (QMA) problem, and 
of random states. The latter example is useful for our 
discussion, since it is known that the average entan¬ 
glement of a random pure state \ifi) oiN qubits is close 
to maximal. Indeed, making use of the von Neumann 


entropy 

'S'(p) =-Tr(plogp), (1) 

it is possible to show mm that the average entropy 
of the reduced state a = Tt:n/2\'4’){'^\ to N/2 qubits 
reads: 

N 

S{a) = — log2 — c -I- (subleading terms), (2) 

where the constant c depends on the choice of the en¬ 
semble used to generate random states mm- For con¬ 
venience we will consider natural logarithm along the 
work. In the case the average is defined with respect 
to the unitarily invariant Haar measure on the space 
of unit vectors of size 2^ the constant reads c = 1/2 
[T] . Although random states are almost maximally en¬ 
tangled, their entropy differs from the maximal value 
by a negative constant. One may then investigate 
when truly maximally entangled states appear, what 
properties they have and what they are useful for. A 
relevant motivation for this study is the role that ab¬ 
solute maximal entanglement may play in the context 
of holography, as we shall discuss later on. 

The aim of the present paper is to extend the rela¬ 
tion of absolutely maximally entangled (AME) states 
to a different branch of mathematics, so-called combi¬ 
natorial designs. On one hand, new families of maxi¬ 
mally entangled multipartite states are presented. On 
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the other hand, we demonstrate a direct link between 
multipartite entanglement and combinatorial struc¬ 
tures, including mutually orthogonal Latin squares 
and Latin cubes, and symmetric sudoku designs. Fur¬ 
thermore, we introduce the concept of multiunitary 
matrices, which exists in power prime dimensions, 
show their usage in constructing AME states and 
present a small catalog of such matrices in small di¬ 
mensions. 


This work is organized as follows. In Section |T^ 
we present some important aspects of entanglement 
theory, AME states and its link with holography. In 
Section m we review the current state of the art of 


AME states for multipartite systems. In Section IV 


we consider the particular case of AME states hav¬ 
ing minimal support and its relation to classical codes 
and orthogonal arrays of index unity. A link between 
AME states, mutually orthogonal Latin squares and 
Latin cubes and hypercubes is explored in Section |Vl 
In Section |VI| we introduce the concept of multiu¬ 
nitary matrices and demonstrate that they are one- 
to-one connected with AME states. Additionally, we 
provide some constructions of multi-unitary matrices 
and its associated AME states. In Section lylTl we re¬ 
sume and discuss the most important results obtained 
in our work and conclude. In Appendix we dis¬ 
cuss the non-existence of AME states of four qubits 
and present the most entangled known states. In Ap¬ 
pendix we detailedly explain the simplest case of 
2-unitary matrices of order D = cP, which are asso¬ 
ciated to AME states of 4 qudits. some examples of 
multiunitary matrices are presented. In Appendix [C] 
we establish a relationship between AME states and 
a special class of Sudoku. Finally, in Appendix [P] we 
present a mini-catalog of multi-unitary matrices ex¬ 
isting in low dimensions. 


II. AME STATES: DEFINITION AND BASIC 
PROPERTIES 

A. Definition 

A lot of attention was recently paid to identify en¬ 
tangled states of A^-party systems, such that tracing 
out arbitrary N — k subsystems, the remaining k sub¬ 
systems are maximally mixed [MHl- Such states are 
often called fc-uniform [HI na and by construction 
the integer number k cannot exceed N/2. In this pa¬ 
per, we shall focus on the extremal case, k = [N/2\ 
(we put the floor function [J to include cases of N even 
and odd) , and analyze properties of states called ab¬ 
solutely maximally entangled (AME) - see [131 El- 
Also, such states were previously known as perfect 
maximally multipartite entangled states |29j . 

The definition of AME states corresponds to those 


quantum states that carry maximum entropy in all 
their bi-partitions. It is a remarkable fact that the 
existence of such states is not at all trivial and deepens 
into several branches of mathematics. Let us be more 
precise and define an AME(A^,(i) state |'0) € H, made 
with N qudits of local dimension d,'H = (C"^)®-^ as a 
state such that all its reduced density matrices in any 
subspace A = , % = A® A, carry maximal 

entropy 


*5'(pa) = y logd VA . (3) 

This is tantamount to asking the reduced density ma¬ 
trices to k qudits to be proportional to the identity 

1 N 

Pk = y/d-" yk< — . (4) 

Let us note the fact that a fc-uniform state is also 
A:'-uniform for any 0 < k' < k. 

There is an obstruction for a state to reach maxi¬ 
mal entanglement in all partitions due to the concept 
of monogamy of entanglement [13 [S]. Every local de¬ 
gree of freedom that tries to get maximally entangled 
with another one is, then, forced to disentangle from 
any third party. Therefore, entanglement can be seen 
as a resource to be shared with other parties. If two 
local degrees of freedom get largely entangled among 
themselves, then they are less able to be entangled 
with the rest of the system. But this rule is not al¬ 
ways fulfilled. There are cases where the values of the 
local dimensions d and the total number of qudits N 
are such that AME states exist. Eor a given N, there 
is always a large enough d for which there exists an 
AME state m- However, the lowest value of d such 
that an AME state exists is not known in general. 

Let us mention that AME states are useful and nec¬ 
essary to accomplish certain classes of multipartite 
protocols. In particular, in Ref. m, it was shown 
that AME states are needed to implement two dif¬ 
ferent categories of protocols. First, they are needed 
to achieve perfect multipartite teleportation. Second, 
they provide the resource needed for quantum secret 
sharing. These connections hint at further relations 
between AME states and different branches of Mathe¬ 
matics. Eor instance, AME states are related to Reed- 
Solomon codes [12!; Also, AME states (and A:-uniform 
states in general) are deeply linked to error correction 
codes m- 

There is yet another surprising connection between 
AME states and holography [T8| [19] . It can be seen 
that AME states provide the basis for a tensor network 
structure that distributes entanglement in a most ef¬ 
ficient and isotropic way. This tensor network can 
be proven to deliver holographic codes, that may be 
useful as quantum memories and as microscopic mod¬ 
els for quantum gravity. A key property for these 





3 


new developments is related to the properties of multi- 


unitarity that will be explored in Sect. VIB 


B. Local Unitary equivalence 

Entanglement is invariant under choices of local ba¬ 
sis. It is then natural to introduce the concept of Local 
Unitary (LU) equivalence among AME states. Two 
quantum states |$) and |5') are called LU-equivalent 
if there exist N local unitary matrices, Ui,... Un such 
that 


1$) = Ui (g) Ua 0 ... 0 . (5) 

If is a maximally entangled state any other state 
LU-equivalent to it is also maximally entangled. In 
this respect, we will define AME(7V, d) as the set of 
all AME states in the Hilbert space 'H{N,d) and de¬ 
note their elements by a Greek letter, e. g. |H 4 _ 3 ) G 
AME(4,3) is an AME of four qutrits. 

The LU transformations introduce equivalence 
classes of states. A question arises naturally about 
which state should be chosen as the representative of 
the class, that will be denoted as canonical form of 
an AME state. It is possible to argue in two differ¬ 
ent directions. On one hand, we may consider that a 
natural representative may carry all the elements of 
the computational basis. It would then be necessary 
to establish theorems and a criterion to fix the coeffi¬ 
cients. On the other hand, an alternative possibility 
is to choose the element of the class with a minimal 
support on the computational basis. Results in both 
directions are presented in Sec. m 

It is not known in general how many different LU 
classes there are in the set AME(iV, d) for every N and 
d. This question can be tackled by the construction 
of LU-invariants. A few examples are at hand for few 
qubits. For three qubits, it is known how to obtain 
a canonical form of any state using LU and that all 
states are classified by 5 invariants [20], only one of 
them is genuinely multipartite, the tangle. For four 
qubits, there are ways to construct a canonical form 
and to find the hyperdeterminant, as well. Yet, it 
is unknown how to proceed to larger local dimensions 
and number of parties. It is arguable that the subset of 
AME states is characterized by several LU-invariants, 
probably related to distinct physical tasks. In such a 
case, there would be different AME states not related 
to each other by LU. 


C. AME and holography 

Quantum holography amounts to the fact that the 
information content of a quantum system is that of 
its boundary. It follows that the information present 


in the system is far less than the maximum allowed. 
Degrees of freedom in the bulk will not carry maximal 
correlations, neither the von Neumann entropy of any 
sub-part of the system will scale as its volume. 

To gain insight in quantum holography, it is natural 
to investigate the bulk/boundary correspondence of 
the operator content of the theory. On the other hand. 
Quantum Information brings a new point of view on 
this issue, since it focuses on the properties of states 
rather than on the dynamics that generates it. In 
this novel context, we may ask what is the structure 
of quantum states that display holographic properties. 
That is, we aim at finding which is the detailed entan¬ 
glement scaffolding that guarantees that information 
flows from the boundary to the bulk of a system in a 
perfect way. 

A concept separate from holography turns out to 
be very useful to address the analysis of holography 
from this new Quantum Information perspective, that 
is Tensor Networks of the kind of matrix product 
states (MRS), projected entangled pair states (PEPS) 
and multi-scale entanglement renormalization ansatz 
(MERA) Indeed, Tensor Networks provide a frame to 
analyze how correlations get distributed in quantum 
states, and thus to understand holography at the level 
of quantum states. Each connection among ancillary 
indices quantifies the amount of entanglement which 
links parts of the system. Holography must necessar¬ 
ily relay on some very peculiar entanglement struc¬ 
ture. 

A first attempt to understand the basic property 
behind holography of quantum states was presented 
in [18]. There, it was proposed to create a quan¬ 
tum state on a triangular lattice based on a tensor 
network that uses as ancillary states an absolutely 
maximally entangled state. To be precise, the state 
\Q) G AME{4:,3) (see Eq.(ll| was defined on tetra¬ 
hedrons, in such a way that the vertices in its basis 
connect the tensor network and the tip of the tetra¬ 
hedron corresponds to a physical index. 

In a subsequent work another construction was 
based on the 5-qubit |T 5 _ 2 ) G AME{5,2) state 
Again, the fact that the internal construction of the 
state is based on isometries is at the heart of the holo¬ 
graphic property. 

There are two obvious observations on the surpris¬ 
ing relation between AME states and holography. The 
first is related to the natural link between AME states 
and error correction codes. It is arguable, then, that 
the essence of holography is error correction, which 
limits the amount of information in the system. The 
second is that the very property responsible for holog¬ 
raphy is multi-unitarity, which is analyzed in depth in 
section [VT] It is further arguable that multi-unitarity 
is the building block of symmetries, since the sense 
of direction is lost and can be defined at will. Those 
ideas deserve a much deeper analysis. 
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D. Related definitions 


1. Maximally Entangled sets 


In the context of convertibility of states via Local 
Operations and Classical Communication (LOCC), 
multipartite entanglement is significantly different 
from the bipartite one. While in the bipartite case, 
there is a single maximally entangled state (up to lo¬ 
cal unitaries) that can be transformed into any other 
state by LOCC (and cannot be obtained from any 
other), in the multipartite scenario this is no longer 
true. In Ref. m, the notion of the Maximally Entan¬ 
gled (ME) set of fV-partite states is introduced as the 
set of states from which any state outside of it can be 
obtained via LOCC from one of the states within the 
set and no state in the set can be obtained from any 
other state via LOCC. Note that this notion of max¬ 
imal entanglement is strictly weaker than the AME, 
in the sense that most of (or all) states in the ME set 
will not be an AME state, but any AME state will be 
in its corresponding ME set. In Ref. [H], the ME set 
is characterized for the case of three and four qubits. 
It is interesting to point out that, unlike the 3-qubit 
case, deterministic LOCC transformations are almost 
never possible among fully entangled 4-partite states. 
As a consequence of this, while the ME set is of mea¬ 
sure zero for 3-qubit states, almost all states are in 
the 4-qubit ME set. This suggests the following pic¬ 
ture; given a fixed local dimension and for an increas¬ 
ing number of parties, the AME states become more 
an more rare at the same time that more and more 
states need to be included in the ME set. In other 
words, while maximally entangled states defined from 
an operational point of view become typical when the 
number of parties increases, AME states are exotic. 

The issue of inter-convertibility between quantum 
states has been addressed under a larger set of op¬ 
erations than the LOCC: the Stochastic Local Oper¬ 
ations and Classical Communication (SLOCC). The 
SLOCC identify states that can be interconverted by 
LOCC in a non-deterministic way, but with a non¬ 
zero probability of success. In this respect, a sys¬ 
tematic classification of multipartite entanglement in 
terms of equivalence classes of states under SLOCC is 
presented in Ref. m- In particular, it is shown that 
the SLOCC equivalence class of multipartite states is 
characterized by ratios of homogeneous polynomials 
that are invariant under local action of the special 
linear group. This work generalizes for an arbitrary 
number of partitions and finite local dimension the 
complete classification made for 4-qubits in Ref. |22) . 


2. Maximally Multipartite Entangled states 

In Ref. [7], Maximally Multipartite Entangled 
(MME) states are introduced as those states that 
maximize the average entanglement (measured in 
terms of purity) where the average is taken over all 
the balanced bipartitions i. e. \A\ = \_N/2\. More 
specifically, the MME states are defined as the mini- 
mizers of the potential of multipartite entanglement 



where = Tr(/93^) is the purity of the partition 
A. Note that the above potential is bounded by 
l/dW2J < T^ME £ 1 and its lower bound is only 
saturated by AME states. 

By minimizing the multipartite entanglement po¬ 
tential, explicit examples of AME states of 5 and 6 
qubits are presented in Ref. [7]. It is remarkable that 
even for a relatively small number of qubits {N > 7), 
such minimization problem has a landscape of the pa¬ 
rameter space with a large number of local minima, 
what implies a very slow convergence. The reason for 
that is frustration. The condition that purity satu¬ 
rates its minimum can be satisfied for some but not 
for all the bipartitions (see [S] for details). In this 
respect, in |24] . the minimization of the multipartite 
entanglement potential is mapped into a classical sta¬ 
tistical mechanics problem. The multipartite entan¬ 
glement potential is seen there as a Hamiltonian which 
is minimized by simulated annealing techniques. 

III. EXAMPLES OF AME STATES 
A. Qubit AMEs 

Let us consider states made out of qubits, that is 
the dimension of the local Hilbert space is d = 2. 
The simplest cases of AME states are any of the Bell 
states. There is a unique partition of two qubits and 
it is possible to entangle both parties maximally. It 
is easy to argue that there is a unique quantity that 
describes the amount of entanglement in the system. 
This can be chosen to be the first eigenvalue of the 
Schmidt decomposition of the state, or some other 
quantity derived from it as the von Neumann entropy. 
All two particle states, whatever their local dimension, 
can be entangled maximally. Those states are of no 
interest for our present discussion which is genuinely 
centered in multipartite entanglement. 

In the case of 3 qubits the well-known GHZ state 
[25] is an AME, but it is known that there is no 4- 
qubit AME j^. The amount of degrees of freedom in 
the definition of the state is insufficient to fulfill all the 
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constraints coming from the requirement of maximum 
entanglement. 

For 5 and 6 qubits, there are AME states. In par¬ 
ticular, a 5-qubit state 1 X 5 ^ 2 ) G AME{5,2) can be 
defined by the coefficients of the superposition of ba¬ 
sis states that form it: 

|T5-2) = 2^i:''r’l>>. (7) 

where we used the usual shorthand notation for the el¬ 
ements in the computational basis and the coefficients 
have the same modulus and signs given by [7] 

= { 1 , 1 , 1 , 1 , 1 ,- 1 ,- 1 , 1 , 1 ,- 1 ,- 1 , 

- 1 ,- 1 , 1 , 1 }. ( 8 ) 

A state AME(5,2) found in [TU] was proven to be use¬ 
ful for a number of multipartite tasks in Ref. [26]. It 
can also be found as the superposition of a perfect er¬ 
ror correcting code as presented in Ref. m and also 
from orthogonal arrays [12]. 

For the sake of completeness, let us also provide an 
absolutely maximally entangled 6-qubit state 

2^ — 1 

1 ^ 6 . 2 ) = (9) 

V 2 2^0 

with 

c(=) = 1,1,1, 

-1,-1,-1,1,1,-1,-1,-1, 

- 1 ,- 1 , 1 ,- 1 ,- 1 , 1 ,- 1 ,- 1 , 

1,1,-1,1,-1,1,-1,-1, 

- 1 , 1 ,- 1 ,- 1 ,- 1 ,- 1 , 1 ,- 1 , 

1 ,- 1 , 1 , 1 ,- 1 ,- 1 , 1 ,- 1 , 

1 ,- 1 ,- 1 ,- 1 , 1 , 1 , 1 ,- 1 , 

1 ,- 1 ,- 1 ,- 1 ,- 1 ,- 1 ,- 1 , 1 }. ( 10 ) 

The case of 7 qubits remains unsolved but numerical 
evidence hints at the impossibility of finding such an 
AME state [28]. Eor 8 qubits or more, there are no 
AME states m- 

B. A central example: AME(4,3) 

The first non-trivial example of AME for larger local 
dimensions corresponds to a state made of 4 qutrits, 
AME(4,3). Its explicit construction is 

1 ^ 4 . 3 ) = ^ l*)b')l* + i)l* + 2j) , (11) 

i,j,=0,1,2 


where all qutrit indices are computed mod(3). It is 
easy to verify that all the reduced density matrices to 
two qutrits are equal to p = ^/g, so that this state 
carries entropy S = 2 log 3 for every one of its bi¬ 
partitions. 

The state 1 114 , 3 ) can be viewed as a map of a two- 
qutrit product basis into a second one. That is 

III4.3) = XI = X (12) 

2 — 0,8 0,8 

where {|ui)} and {|ui)} are product basis for two 
qutrits, and |ui) = Uij\uj). The matrix Uij is not 
only unitary (as it must as a consequence of multiu- 
nitarity) but also a bijective map between the sets 
of words of length 2 over an alphabet Z 3 = {0,1,2}. 
In other words, the entries of Uij are 0 with a single 
1 per row/column. This property remains the same 
whatever partition is analysed, though the unitary will 
vary. 

A second feature of the state 1 114 , 3 ) is that the Ham¬ 
ming distance between any pair of elements in the 
state is three, Dh = 3. As all the sequences in 1 114 , 3 ) 
differ in 3 elements, any single qutrit error can be cor¬ 
rected. Hence, maximal entanglement is related to 
error correction codes. Both of these properties are 
related to the fact that 1 114 , 3 ) is an AME of minimal 
support, a concept we are going to explore in Sect. 

Cyj 

It is also natural to expect AME states to accom¬ 
modate easily to some magic-square like relations. A 
simple example goes as follows. Write the coefficients 
of 1114 , 3 ) as a matrix of row i and column i, giving the 
composed value a of the remaining two qutrits from 
0 to 8, that is a = 3(i -I- j)mod(3) + {i + 2j)mod(3). 
The square reads 

0 5 7 
4 6 2 
8 1 3 

where all rows and columns add up to 12. The same 
properties are maintain if we interchange the indices 
in the state. These kind of combinatorial designs are 
going to be explored in sect. |Vj 

From an experimental point of view, 1 114 , 3 ) can be 
created defining a quantum circuit that generates it. 
Such a circuit makes use of the following two gates: 

Fourier F 3 IO) = ^ (|0) -b |1) -b |2)), 

C 3 - adder Uc-adder\i)\j) = |i)|(i-b j)mod3). 

(13) 

This gate C 3 generalizes the CNOT gate for qubits 
and it is represented in the circuit using the usual 
symbol for CNOT with the subscript 3. The |ll 4 , 3 ) 
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the case of a certain bipartition j of the set 5'". The 
real coefficients are free as long as they satisfy 
Eq. (15) and the normalization condition. 

With this expressions all AMEs of a small number of 
qubits are classified. It is also worth noting that AME 
states with maximal support and uniform amplitudes 
Tfc = \/k are possible to define (by just setting 

all Cj to 0 in Eq. ®) m- 


FIG. 1: Quantum circuit required to generate the state 
AME fl 4,3 (4 qutrits) based on the Fourier gate fa and 
control-adders gate Ca — adder . 

state can be constructed as a sequence of a Fourier F 3 
and Ca-adders following the circuit depicted in Fig[l] 
acting on the initial state |0000) of 4 qutrits. 

C. General expression for AME states 

The states AME(5,2) and AME(6,2) have the maxi¬ 
mal number of terms whereas AME(4,3) has the mini¬ 
mal possible number. Therefore the following question 
arises: is there a general expression for AME(7V, d) 
having the maximal number of terms? Facchi par¬ 
tially solved this issue for qubits in [^. It can be 
shown that a state expressed as 

2"-l 

w = E Zk\k), Zk=rk^k, (14) 

k—0 

where the \k) span the whole computational basis and 
Tfc and are respectively the modulus and the phase 
of the complex coefficients Zk, is an AME if it satisfies 
the following equations: 

G 0 mG' 0 m^l 0 m^l' © m 0 ; (15) 

m 

where I, V and m stand for both parts of a certain 
balanced bipartition (of [A^/2J particles). 

The general form of the squared modulus of the 
coefficients will be: 

ri=2-«+ •£ E4"’ 

^<n<N l<h<n 

where S'" stands for the set of bipartitions of the sys¬ 
tem into groups of n and N — n particles and j is an 
index for the bipartitions of this set. h is an index 
for particles contained in the bipartition and 
stands for the value (0 or 1 ) of particle h in ket k, in 


IV. AME STATES OF MINIMAL SUPPORT 
AND CLASSICAL CODES 

In Ref. [30], a subclass of AME(7V, d) states is shown 
to be constructed by means of classical maximum- 
distance separable (MDS) codes. In this section we 
show that such a subclass corresponds to the set of 
AME states of minimal support and exploit these 
ideas to get conditions for their existence. 


A. Support of AME states 

From the explicit examples we have presented in 
Sec. |III[ AME states appear to need a different 
numbers of elements to be written. For instance, 
AME(4,3) is made of the superposition of 9 states, all 
weighted with the same coefficient. Yet, AME(6,2) 
as written using the 64 basis states with coefficients 
either 1 or - 1 . 

Let us define the support of a state |^) as the num¬ 
ber of non-zero coefficients when Itjj) is written in the 
computational basis. The support of a class is defined 
as the support of the state inside the class with mini¬ 
mal support. Note that the support of a class defines 
in turn another equivalence class. Two states are sup¬ 
port equivalent if they belong to a LU classes with 
equal support. 

In this sense, it is interesting to point out that the 
state AME(6,2), defined in Eq. Q with the maximal 
support of 2® states, is LU equivalent to a state of sup¬ 
port 16 by applying some Hadamard gates on its basis. 
It can be proven that 8 is the minimal support that 
such state could have but it is not attainable. Also in 
[To] a state AME(5,2) is built using 8 elements, while 
the theoretical minimum would be 4. We may won¬ 
der why the naive minimum possible number of 2 L^/^J 
elements is not always attained. The answer to this 
question is answered in Ref. [12] . where it is proven a 
one-to-one relationship between /c-uniform states hav¬ 
ing minimal support and a kind of combinatorial ar¬ 
rangements known as orthogonal arrays. Therefore, 
the non-existence of such states having minimal sup¬ 
port is due to the non-existence of some classes of 
orthogonal arrays (those having strength one). In the 
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following subsections we study in detail AME states 
having minimal support. 

B. Equivalence between AME states of minimal 
snpport and MDS codes 

A maximally entangled state in AME( A, d) belongs 
to the class of minimal support iff it is LU-equivalent 
to a state |t[') with support i. e. 

dLiV/2J 

= (17) 

k=l 

where Xk G are words of length N over the alpha¬ 
bet Zd = { 0 ,..., d — 1 }, and > 0 and 9k € [ 0 , 27r) 
are their modulus and phases respectively. 

Given a bipartition A = {oi,..., a|^„j }, it will 
be useful to introduce the subword Xk[-A\ of the 
word Xk for the partition A as the concatenation 
of the ai-th, a 2 -th, ... a„-th letters of Xk, that is, 
Xk{A] = a;fc[ai]a;fc[a 2 ]... Xk[an]- Let us also denote by 
= {xk,k = 1 ,..., } the set of words which 

14') has support on, and by = {xk[A\,k = 

1 ,..., dL^/^J } the set of all subwords Xk[A\ corre¬ 
sponding to the bipartition A. With this notation, 
the reduced density matrix of a partition A can be 
written as 

PA^^{xk'\A]\\^)m\xk'[A]), ( 18 ) 

k' 

where 


differ, e. g. d//(00010,10000) = 2. To see this, note 
that otherwise there would exist a balanced biparti¬ 
tion A' for which Xi \A'] = Xj [A'] for i ^ j and conse¬ 
quently the set \A'] would not contain all the pos¬ 
sible words of length [1V/2J, that is, C . 

A set of M words of length N over an alphabet 
of size d that differ pairwise by at least a Hamming 
distance 5 is called a classical code. What we have 
shown above is that the existence of AME( A, d) states 
of minimal support imply the existence of classical 
codes of M = |A,i,| = words with Hamming 

distance 6 = [N/2\ + 1. The codes produced by AME 
states are special in the sense that they saturate the 
Singleton bound EB, 

M < d^-^+^. (21) 

This type of codes that saturate the Singleton bound 
are called maximum-distance separable (MDS) codes. 

The converse statement is also true. That is, the 
existence of a MDS code also implies the existence of 
an AME state with minimal support [5D]. The argu¬ 
ment is the following: a code of M = words 

of length N and Hamming distance 5 = [1V/2J + 1 
has all its subwords associated to any balanced bipar¬ 
tition A of size \A\ = [A^/2J different, which implies 
the condition 3 above. Thus, AME states of minimal 
support are equivalent to classical MDS codes. 

This equivalence can be exploited to see that a nec¬ 
essary condition for AME(A, d) states with minimal 
support (and equivalently of MDS codes) to exist is 
that the local dimension d and the number of parties 
N fulfill 


= ThW H rke^^'‘\xk[A])\xk[A]). (19) 

k=l 

In order for |'I') to be an AME, the reduced density 
matrix of any bipartition A = {ai,..., a|^jv/ 2 j } needs 
to be the completely mixed state. It is easy to see 
that this has the following implications: 

1. The modulus Vk = 1 for all k. 

2. The phases 9k are arbitrary. 

3. Eor any balanced bipartition A, with |A| = 
[IV/2J, two words Xi , Xj G Xqi have subwords 
Xi[A\ = Xj[A\ if and only \i i = j. Equivalently, 

Condition 3 implies that any pair of different words 
Xi , Xj G have Hamming distance 

dff(x„x,) > LA/2J+1, (20) 

where the Hamming distance between two codewords 
is defined as the number of positions in which they 


d > [7V/2J -f 1. (22) 

We can prove it as follows: let us try to construct 
an AME state by building an MDS code. Due to the 
relabeling freedom, the first d -I- 1 words of the code 
can be chosen as 


d+ 1 
code 
words 


0 

0 


0 

. 0 


. 00 0 . 0 , 

. 01 1 . 1 , 

...O(d-l) (d-l)...(d-l), 

.10 Xd[[N/2\+l]...Xd[N], 


where the letters Xd[i]: for [1V/2J -1-1 < i < A, are still 
unknown and every word is written in two subwords 
of lengths [iV/2j and [(A-|-1)/2J respectively. 

Note that: (i) none of the unknown letters Xd[i] 
can be 0 in order for the word Xd to have Hamming 
distance L 1 V/ 2 J -|- 1 with the first word 0 ... 00 ... 0 . 
{ii) None of the unknown letters Xd[i] can be repeated 
in order for Xd to have Hamming distance [1V/2J -|- 1 
with the other d — 1 words. Therefore, if [{N + 1)/2J 
variables, Xd[i], must take [(iV-|-l)/2j different values. 
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and all of them must be different from 0 , it is necessary 
to extend the alphabet, forcing that d > \N/2\ + 1 . 

Interestingly, Eq.(22) forbids the existence of 
AME(N,2) states having minimal support for N > 
3. However, this does not represent a proof that 
AME(4,2) do not exist. Also, this inequality is sat¬ 
urated for the states AME(4,3) and AME(6,4). The 
existence of the cases AME(8,5) and AME( 8 , 6 ) is still 
open whereas the states AME(8,7) and AME( 8 , 8 ) are 
known. 


C. A less trivial example: AME(6,4) 


An application of the above connection between 
AME states and MDS codes is the construction of 
AME states by exploring the set of all words and se¬ 
lecting those which differ in at least [A^/2J -|- 1 ele¬ 
ments. For the case of the AME(6,4), such a search 
gives a state with a equal superposition of the follow¬ 
ing entries: 

{ 000000,001111,002222,003333,010123,011032, 
012301,013210,020231,021320,022013,023102, 
030312,031203,032130,033021,100132,101023, 
102310,103201,110011,111100,112233,113322, 
120303,121212,122121,123030,130220,131331, 
132002,133113,200213,201302,202031,203120, 
210330, 211221, 212112, 213003, 220022, 221133, 
222200, 223311, 230101, 231010, 232323, 233232, 
300321,301230,302103,303012,310202,311313, 
312020,313131,320110,321001,322332,323223, 
330033,331122,332211,333300} . (23) 


Such an AME state carries the minimum possible sup¬ 
port. 


D. Construction of AMEs with minimal support 

Finding AME states by exploring the set of all 
words is highly inefficient and becomes in practice un¬ 
feasible from a relatively small number of parties. In 
this context, the Reed-Solomon codes m can be a 
useful tool to produce systematic construction of MDS 
codes and equivalently AME states. 

Let us review here the particular case of d prime 
and N = d + 1. Let us refer to the elements of the 
superposition in the quantum states as words Xi and 
the word of a half-partition as ui. The code words are 
obtained using the action of a generator G, Xi = Ui- G. 
The problem is then reduced to fixing G. It can be 


shown that a family of valid generators is given by 

/ 1 1 ... 1 0 \ 

^ _ 9o 9i ■■ ■ 9d-i 0 - 


V 5o 5? ■ ■ ■ 9Li 1 

where d S Prime, IV = d -|- 1, and k = N/2. 

The case of AME(4,3) can be re-obtained using 

«=((!)!) 

Another concrete example corresponds to go = !> 
gi = 1,..., go = 6 that corresponds to AME(8,7), with 
N = 8 d = 7-dits. Then a total of 7'^ codewords are 
obtained that differ by a minimum Hamming distance 
dH = 5 . 

As the above construction can only be accomplished 
for d being prime and N = d+1, it is interesting to ad¬ 
dress the question of whether given some AME (TV, d) 
it is possible to construct another AME(A^',d'). In 
this context, the following lemma can be useful which 
allows to construct an AME(iV',d) for any N' < N: 
if there exists an AME state with minimal support in 
7i{N,d) with N being even, then there exist other 
AME states with minimal support in the Hilbert 
spaces HiN'^d) for any N' < N. 

To prove this last statement, we will consider sepa¬ 
rately the transitions N ^ N — 1 and — 1 —>■ N — 2. 

Transition N —>■ iV — 1: The existence of an AME 
state with minimal support implies the existence of a 
code of d^/^ words of length N with Hamming dis¬ 
tance dn > A^/2 -I- 1. Let us order the words in the 
code in increasing order and take the subset of the 
first words which start with 0. Note that by 

suppressing such a 0, we get a code of words 

of length — 1 with Hammings distance N/2 + 1, 
forming an AME. 

Transition A^ — 1 —)■ N — 2: From the previous steps, 
we are left with a code of words of length 

A^ — 1 and Hamming distance N/2 + 1. Note that 
by suppressing an arbitrary letter from all the words 
of the code, one is left with a set of words of 

length N — 2 with Hamming distance A/2-1-1 — 1 = 
{N — 2)/2 -I- 1, which is a MDS code. 

By iterating the previous procedure, we obtain 
MDS codes (and AME states of minimal support) 
for any 2 < N' < N. Equivalently, MDS codes 
are constructed from considering orthogonal arrays of 
strength one [T^ . 


E. Non-minimal support AMEs and perfect 
quantum error correcting codes 

AME states are deeply related to classical error cor¬ 
rection codes and compression US HZl El]. This is 
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somewhat intuitive since maximal entropy is related 
to maximally mixed subsets. The measure of any local 
degree of freedom delivers an output which is com¬ 
pletely random. This is, in turn, the basic element 
to correct errors. Hence, a relation between the el¬ 
ements superposed to form an AME state and error 
corrections codes is expected. 

Le us illustrate the connection of an AME(5,2) state 
with the well-known 5-qubit code m- It is easy to 
see that by applying some Hadamard gates on local 
qubits, the AME(5,2) state, defined through the 32 
coefficient given in Eq. is LU-equivalent to a state 
with fewer non-zero coefficients. Actually, a represen¬ 
tative of the same AME(5,2) class is found to have 
only 8 coefficients. That state corresponds to a super¬ 
position of the two logical states in the error correcting 
codes found in m 

|f^5.2) = ^(|0U + |lU), (26) 

where the logical qubits are defined as 

|0)l = ^ (100000) -H 100011) -H 101100) - 101111)) 

|l)i = ^ (IllOlO)-H IllOOl)-H IIOIIO) - 110101)) 

Note the fact that the coefficients carry both plus and 
minus signs as is the case in the non-minimal support 
AMEs. 


V. COMBINATORIAL DESIGNS AND 
fe-UNIFORMITY 

Combinatorial designs are arrangements of elements 
satisfying some specific properties so called balanced 

[33] . Such elements are restricted to a finite set, typi¬ 
cally considered as subsets of integer numbers. Some 
remarkable examples are block designs, t-designs, or¬ 
thogonal Latin squares and orthogonal arrays (see 

[34] and references therein). Combinatorial designs 
have important applications in quantum physics |35| . 
Indeed, a connection between genuinely multipartite 
maximally entangled states and orthogonal arrays has 


been recently found |12] . Furthermore, they are a 
fundamental tool in optimal design of experiments 
[361 El]. The existence of some combinatorial de¬ 
signs can be extremely difficult to prove. For example 
the existence of Hadamard matrices in every dimen¬ 
sion multiple of four (i.e., the Hadamard conjecture) is 
open since 1893 [38] and it represents one of the most 
important open problems in Combinatorics. 

A. Relation to mutually orthogonal Latin 
squares 


Let us consider the explicit expression 
AME(4,3) state defined in Eq.(ll): 


of the 


|H4_3) = i (10000) -f 10112) -h 10221) 

o 

-H|1011) -h |1120) -f |1202) 

-f|2022)-h |2101)-f |2210)). (27) 


Here, the third and fourth symbols appearing in every 
term of the state can be arranged into 2 Greco-Latin 
squares of size three: 


Aa Bj Cp 
Bfj Ca Ay 
Cy A/3 Ba 


KJi, go 

KO QlH A* , (28) 

g* AO 


where every symbol of the sets {A, a}, {B,/3} and 
{C, y}, is associated to 0,1 and 2, respectively. Hence 
a pair of symbols may represent a card of a given rank 
and suit. Note that the first two digits in Eq.(27) may 
be interpreted as addresses determining the position of 
a symbol in the square. Furthermore, by considering 
the following four Mutually Orthogonal Latin Squares 
(MOLS) [39] of size five 


0000 

4321 

3142 

2413 

1234 

1111 

0432 

4203 

3024 

2340 

2222 

1043 

0314 

4130 

3401 , 

3333 

2104 

1420 

0241 

4012 

4444 

3210 

2031 

1302 

0123 


(29) 


we define a 2-uniform state of 6 subsystems with five 
levels each: 


|$5) =i (1000000)-h 1104321)-H 1203142)-h 1302413)-H 1401234)-H 
1011111) -f 1110432) -h 1214203) -h |313024) -h |412340) -f 

1022222) -h 1121043) -h |220314) -h |324130) -f |423401) -f 

1033333) -h 1132104) -h |231420) -h |330241) -h |434012) -h 

1044444) -h 1143210) -h |242031) -h |341302) -h |440123)). 


( 30 ) 
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A state locally equivalent to j^s) has been previously 
found in m but its connection to MOLS is firstly 
given here. By considering the standard construction 
of maximal sets of d — 1 MOLS of prime size d we can 
generalize the above construction for quantum states 
of a prime number of levels d and TV = d + 1 parties 
as follows 


^ d—1 d—1 

i,j—0 m—1 


d = 10 expression (33) describes a 2 -uniform state of 
4 subsystems with 10 levels each. 

In general, the problem of constructing TV—2 MOLS 
of size d is equivalent to construct a 2 -uniform state 
of TV qudits of d levels having d^ positive terms. Note 
that d^ is the minimal number of terms that can have 
a 2 -uniform state of qudits of d level systems. 


B. AME states and hypercubes 


It is well known that a maximal set of d — 1 MOLS 
of size d exist for every prime power d = p"* uni- 
This means that the above general expression can be 
extended to the case of prime power level systems. For 
instance, the maximal set of 3 MOLS of order d = 4 
can be represented by a color figure 


A4 

A'O 

go 


TO 

g<ik 

a:* 

AO 

g* 

JO 

AO 

Kilt 


AUt 


go 


This design determines an AME(5,4) given by 

105 , 4 ) = 100000) + 110312) -k 120231) -f |30123) -f 
101111) -f 111203) -f 121320) -f 131032) -k 
102222) -h 112130) -f |22013) -f |32301) -f 
103333) -f 113021) -k |23102) -k |33210), 

where every symbol of the sets {A, 4|k, blue}, 
orange}, {Q, Jit, green} and {AT, O, red}, is as¬ 
sociated to 0,1, 2 and 3, respectively, while the first 
two digits of every term label the position of a symbol 
in the pattern. In the above expression a normaliza¬ 
tion factor is required. Note that this state, or a state 
equivalent with respect to local unitary transforma¬ 
tions, arises from the Red-Solomon code of length five 

m- 

Furthermore, the construction can be extended to 
any dimension d in the following way: 

d-l N(d) 

|^d-Hl,d) = )) X! I*’-)) 0 \^rn[i,j]), (33) 

i,j—0 m—1 

where TV(d) denotes the maximal number of MOLS 
of size d. Here \m\ij] denotes the entries of the m-th 
Latin square, so the above expression can be consid¬ 
ered as a direct generalization of Eq.(2). It is worth 
to add that for dimensions d > 12 not equal to a 
prime power number only lower bounds for the func¬ 
tion TV(d) are known [SU]- The problem is solved only 
for smaller dimensions, as Af (6) = 1 in agreement with 
unsolvability of the famous Euler problem of 36 offi¬ 
cers, while TV(10) = 2 - see [40] where explicit form 
of a pair of MOLS of size ten is derived. Thus for 


In the previous section we considered maximal sets 
of MOLS to construct 2-uniform states of qudits. 
However, this construction is not useful to find AME 
states for d > 4. The aim of this section is to 
consider combinatorial arrangements for constructing 
AME states in such cases. The main result is in¬ 
spired in a generalization of the AME(4,3) state jO) 
given in Eq.(27) and the AME state of 6 ququarts 
presented in Eq. (23). In [T^ it was shown that this 


000 

111 

222 

333 

123 

032 

301 

210 

231 

320 

013 

102 

312 

203 

130 

021 


132 

023 

310 

201 

011 

100 

233 

322 

303 

212 

121 

030 

220 

331 

002 

113 


213 302 031 120 
330 221 112 003 
022 133 200 311 
101 010 323 232 


321 230 103 012 
-202 313 020 131 
110 001 332 223 
033 122 211 300 


FIG. 2: Three mutually orthogonal Latin cubes of dimen¬ 
sion 3 and size 4. This arrangement allows us to generate 
a state AME(6,4) of 6 ququarts. Each of the 12 planes 
(4 horizontal, 4 vertical and 4 oblique) contains a set of 3 
MOLS of size 4. 


state can be derived from the irredundant orthogonal 
array IrOA(64,6,4, 3). Eurthermore, this orthogonal 
array can be interpreted as a set of three mutually 
orthogonal Latin cubes (see Eigj^. Also note that 
the AME(4,3) state jO) arises from an IrOA(9,4,3,2) 
(See Eq.(Bl) in (Hj). Thus, if k mutually orthog¬ 
onal hypercubes of dimension k having fc -|- 1 sym¬ 
bols exist they are one to one connected with an 
IrOA((fc -|- lY,2k,k + l,fc) and, therefore, it would 
produce an AME(2k,k-|-l) state. This family of states 
would saturate the bound 


d> y + 1 , 


(34) 


already defined in Eq.(22 1 Note that for fc = 1 one ob¬ 


tains the standard Bell state, |4)i) = (|01) -|- |10))/-\/2 
(1 Latin square of size 2 with 2 symbols). Taking 
fc = 2 and using 2 MOLS of size 3 and 3 symbols we 
arrive to the AME(4,3) state of |H) Eq. 0 which, 
from this point of view, can be considered as a gen¬ 
eralization of the Bell state. Furthermore, the state 
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( [^ corresponding to fc = 3 also belongs to this fam¬ 
ily and it is associated to three mutually orthogonal 
Latin cubes of dimension 3 and size 4. It is interesting 
to check, whether there exist other states with fc > 5 
belonging to this family. 


VI. AME STATES AND MULTI-UNITARY 
MATRICES 

A. Unitary matrices and bipartite systems 

Let us illustrate the connection between unitary 
matrices and AME states for the simplest case of 2 
qubits. Let us assume that the state of the system is 
given by 

10) = ^(C/o,o|00) + t7o.i|01) + 17i,o|10) + Ui,i|ll)), 

(35) 

where pA = \UU\ pb = and T denotes 

transposition. Any unitary matrix U of size 2 repre¬ 
sents a Bell-like state. Furthermore, the Pauli set of 
four unitary matrices lA = {D, (Ta,, cry, ctz}, orthogonal 
in the sense of the Hilbert-Schmidt product, defines 
the maximally entangled Bell basis in 7^2 ® 772 ■ 


B. Multi-unitary for AME(4,3) 

We shall consider again the AME(4,3) state of four 
qutrits 1124 , 3 ), and represent its coefficients by a four- 
index tensor 

1 ^ 24 , 3 ) = E ^ijkl \ijkl). (36) 

i,j,k, 1 = 0 , 1,2 

where the entries of the tensor t can be expressed as 
product of the Kronecker delta functions 

^ijkl — '^^k,i+j^l,i+2j • (^'^) 

Here, the addition operations are modulo 3. As dis¬ 
cussed in the previous section, all non-zero coefficients 
are equal. The tensor tiju consists of 3^ = 81 ele¬ 
ments which can be reshaped to form a square matrix 
of order 9. Note that there exist altogether ( 2 ) = 6 
different ways of choosing a bi-partition of the indices 
and forming a matrix L0,i/. That is, 

( (i 4- 27 , k 21 ), {k 4 - 2 /, i 4 - 2 j) 

{pjiy) = < (z 4-2/c,7 4-21), (j 4-21,7 4-2fc) . (38) 

(1 4- 21, j 4- (7 4- 21 l, i 4- 21) 

The non-trivial property of an AME(4,3) tensor is 
that these six matrices are unitary. As transposition 
of a unitary matrix remains unitary, it is sufficient, 


in this case, to check unitarity for the three cases ap¬ 
pearing in the first column of the right side of Eq. (38). 


That is, taking combined indices in the original tensor. 


^ijkl — U, 


( 1 ) 

(b)(feO 


= 


( 2 ) 

iik){jl) 


= u, 


( 3 ) 

iij){kl)- 


(39) 


absolute maximal entanglement is achieved if the ma¬ 
trices and correspond to different 

changes of bases, that is unitary matrices. 

We refer to such particular kind of unitary matrices 
as multi-unitary. See Appendix for further expla¬ 
nations about reorderings of elements in matrices of 
square size D > A. 


C. General multi-unitarity 

Let us consider a more general case of pure states 
of N subsystems with d levels each. That is, 

d-l 

k) = E (40) 

So,..-,SiV-l=0 

Let us assume here that the number of subsystems 
is even, N = 2k, so there exist M = possible 
splittings of the system into two parts of the same 
size. 

A necessary condition for |0) to be an AME state 
is that the tensor t with 2k indices, reshaped into a 
square matrix of size forms a unitary matrix U. 
This is so, as the reduction associated to the first k 
qudits, given by pk = UU\ should be proportional to 
the identity. To arrive at an AME(2k,k) state, simi¬ 
lar conditions have to hold for all M different square 
matrices obtained from the tensor t by all possible 
ways of reshaping its entries into a square matrix. 
This observation provides a clear motivation to in¬ 
troduce the notion of multi-unitarity: A square ma¬ 
trix A of order (k > 2), acting on a composed 
Hilbert space 77f ^ and represented in a product basis 
by Am,...,nk := (ni,..., nfc|A|z/i,..., z^fe) is Ic-unitary 

if it is unitary for M = reorderings of its entries 
corresponding to all possible choices of k indices out 
of 2k. 

In this way, we can establish the following one-to- 
one connections: 

AME{2, d) = unitary of order d, 

for bipartite systems of 2 qudits having d levels each 
and, in general: 

AME{2k, d) = A:-unitary of order d^, 

for multipartite systems oi N = 2k qudits having d 
levels each. By construction, 1-unitarity reduces to 
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standard unitarity. Any fc-unitary matrix with k > \ 
is called multi-unitary. It is well-known that unitarity 
of matrices is invariant under multiplication. Multi- 
unitarity imposes more restrictions on a given matrix 
[7 than unitarity. Therefore, the product of two multi¬ 
unitary matrices in general is not multi-unitary. For 
instance, the matrix Og (see Eq.(|4l]) below) is hermi- 
tian and 3-unitary, but Og = D is only 1-unitary, as 
it represents a 6-qubit quantum state equivalent to 
GHZ. 


I/T 2 ^ PerTO(0,5,7,l,3,8,2,4,6) and = 

Ferm(0, 2,1,4, 3, 5,8, 7,6), where T 2 and R mean par¬ 
tial transposition and reshuffling (see appendix for 
further details). Here, Perm denotes a permutation 
matrix. A mini-catalogue of all the multi-unitary ma¬ 
trices defined in this work can be found in Appendix 

m 


D. AME and Hadamard matrices 


Similarly, the case of the AME(4,3) state |H) re¬ 
duces to analyzing the properties of the tensor tijki in 
Eq. (36) and verify the multi-unitarity of U. Indeed, 


for this state we have U = Perm{0, 5, 7,4, 6, 2,8,1, 3), 


For six qubits, the AME(6,2) state IS 6 . 2 ) of Eq.(10l 
having a maximal number of terms arises from graph 
states [12]. Let us write it explicitly 


- 6 . 2 / 


I ( -1000000) - 1000001) - 1000010) -I-1000011) - 1000100) -f 1000101) -1-1000110) -b 1000111) 

-1001000) - 1001001) - 1001010) -f 1001011) -b 1001100) - 1001101) - 1001110) - 1001111) 

-1010000) - 1010001) -b 1010010) - 1010011) - 1010100) -b 1010101) - 1010110) - 1010111) 

-blOllOOO) -b 1011001) - 1011010) -b 1011011) - 1011100) -b 1011101) - 1011110) - 1011111) 

-IIOOOOO) -b llOOOOl) - 1100010) - 1100011) - 1100100) - 1100101) -b 1100110) - 1100111) 

-b|101000) - llOlOOl) -b 1101010) -b 1101011) - 1101100) - 1101101) -b 1101110) - 1101111) 

-b|110000) - 1110001) - 1110010) - 1110011) -b 1110100) -b 1110101) -b 1110110) - 1110111) 

-b|111000) - 1111001) - 1111010) - 1111011) - 1111100) - 1111101) - 1111110) -b 1111111)) 


This state leads to the following orthogonal matrix of 
order D = 2^ = 8 which is 3-unit ary: 


Og 


1 

VI 


/ -1 -1 -1 1 

-1 -1 -1 1 

-1 -1 1 -1 

11-11 
-1 1 -1 -1 

1-111 
1 -1 -1 -1 
V 1 - 1-1 -1 


-1 1 1 1 \ 

1 -1 -1 -1 

-1 1 -1 -1 

-1 1 -1 -1 

-1 -1 1 -1 

-1 -1 1 -1 

111-1 
-1 - 1-1 1 / 


(41) 

Note that the entries of |S 6 , 2 ) are given by the con¬ 
catenation of the rows of Og, up to normalization. 
This matrix is symmetric and equivalent up to enphas- 
ing and permutations |43j to the symmetric Hadamard 
matrix Hg = 

Note that Og is 3-unitary but is not, so permu¬ 
tation or enphasing of a unitary matrix can spoil its 
multi-unitarity. Moreover, from the concatenation of 
the rows of 77®^ we only generate a 1-uniform state, 
which means that is only 1-unitary (i.e., unitary). 

We conjecture that for any AME state one can 
choose suitable local unitary operations such that it 
is related to a multi-unitary complex Hadamard ma¬ 


trix. It represents a maximally entangled state with 
the maximal number of terms having all entries of the 
same amplitude. For example, the AME states aris¬ 
ing from coding theory (nidT] and graph states |42| 
are of this form. We recall that fc-uniform states of 
N qudits with d levels having the minimum number 
of terms (d^) are closely related to linear MDS codes 
(See Section 4.3 of [34]) and also to orthogonal arrays 
of index unity |12] . This reasoning implies that any 
pure state of n subsystems with d levels each having 
at least one reduction to n/2 qudits maximally mixed 
can have all its entries of the form if d is even. 


E. Further constructions of AME states 


Note that the state AME(4,3) |H 4 _ 3 ) is not equiv¬ 
alent with respect to local unitaries to a real state 
having all its entries of the form ±3“^. This is a 
consequence of the fact that a real Hadamard matrix 
of size 9 does not exist. However, the state |H 4 _ 3 ) is 
equivalent under local unitary operations to the state 
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(See Eq.(3) in [H]): 


Pm) = 5 E 


j{i-k)+l(i+k) 


(42) 




where uj = This state is associated to the fol¬ 

lowing 2-unitary complex Hadamard matrix: 


/ 


^-=3 


1 
1 
1 
1 

w 

9 

W 

1 

9 

W 

\ w 


1 

1 

1 

w 

9 

W 

1 

2 

W 
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(43) 

Interestingly, every integer power (Up)"^ is a complex 
Hadamard matrix for m ^ 4 (Mod 4) and (C/p)® = D. 
Moreover, Up is equivalent to the tensor product of 
Fourier matrices F 3 ® Fg, that is 


Up = VF3 (g) FlVV, 


(44) 


where V = Diag(l, 1,1, l,w,a;^, l,a;^,a;) is a diagonal 
unitary matrix, while F is a permutation matrix which 
changes the order of the columns from {1,... ,9} to 
{1,4,7,2,5, 8 ,3, 6 ,9}. 

In order to construct a 2-unitary matrix one has to 
take a unitary U such that its partially transpose C/^^ 
and the reshuffled matrix U^ are unitary - see Ap¬ 
pendix 1^ In the case of a matrix U of size D — 3^ 
this implies that the set of nine 3x3 unitary matrices 
appearing in the 3x3 blocks of Eq. (431 define an or¬ 


thogonal basis for the Hilbert-Schmidt product. It is 
thus possible to obtain AME states by considering or¬ 
thogonal bases of unitary operators. For instance, one 
can construct the |H) state from the following matrix: 
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(45) 


We applied here the orthogonal basis defined by the 
displacement operators of size d = 3: 


(46) 

where p = (pi,P 2 ) € 'p = 

X\k) = |fc -I- 1) and Z\k) = uj^\k). These operators 


define the discrete Weyl-Heisenberg group. This ap¬ 
proach can be easily generalized to any prime d > 2. 
Indeed, for d prime every reordering of indices lead us 
to the same matrix up to permutation of columns and 
rows, and therefore it remains unitary. This shows 
a construction of AME(4,d) working for any prime 
number of levels d. Moreover, it is likely that this 
construction can be generalized for prime powers d by 
considering the theory of Galois fields. Observe that 
the above construction is essentially different from the 
construction of AME states used in coding theory. In¬ 
deed, the tensor products of N displacement operators 
bases of size d, i.e. the set {Dpi pi ® ® DpN pN}, 

produce codes and states AME{N,d) [TT] . 


VII. CONCLUSION 


In this work we have analyzed some new proper¬ 
ties of Absolutely Maximally Entangled (AME) states 
in multipartite systems. First of all, we have re¬ 
viewed and extended several ways of constructing 
them. Then, we have explored their relation to 
the field of combinatorial designs. For instance, a 
state AME(4,3) consisting of four maximally entan¬ 
gled qutrits is linked to the set of two mutually orthog¬ 
onal Latin squares of order 3, while a state AME(6,4) 
made of 6 ququarts is related to the set of three mu¬ 
tually orthogonal Latin cubes of order 4. 

A profound relation between AME states and ten¬ 
sors that display the property of multi-unitarity has 
been found: AME states made out of an even num¬ 
ber N of degrees of freedom are equivalent to multi¬ 
unitary matrices (i.e., matrices being unitary after 
M = {j^/ 2 ) rearrangements of its entries). This re¬ 
markable property may be at the core of the use of 
AME states in holography nmn]. 

Furthermore, making use of a state AME{2k,d) con¬ 
sisting oi N = 2k parties, one can construct a fc — 1- 
uniform state state of 2fc — I parties by removing a 
single subsystem. 

We have proven the existence of states AME(4,(i) 
for every prime d > 2. Note that this also provides the 
existence of states AME(3,d) for every prime d > 2. 

AME states remain mainly unexplored. Let us 
bring a number of open problems that deserve to be 
solved. 

• 7 qubits. 

Is there an AME state of 7 qubits? 

• Classification of AME and LU invariants. 

It is unclear whether there is a clear cut classi¬ 
fication of AME, which is related to LU invari¬ 
ants. The fact that some AME states carry dif¬ 
ferent minimal supports or that some AME are 
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right away related to Reed-Solomon codes hints 
a some unknown structure among AME states. 

• Non minimal support. 

Examples of AME states with non-minimal sup¬ 
port are only explicitly known for 5 and 6 qubits. 
It would be natural to find examples for higher 
local dimensions. 

• Computation of invariants. 

LU invariants grow exponentially with the size 
of the number of parties. An example of them 
is the hyperdeterminant, which has only been 
computed up to 4 qubits. There are no compu¬ 
tations of hyperdeterminants of 4 qutrits. Do 
AME states carry maximum values for some LU 
invariant? Erom the general theory of hyperde¬ 
terminants, we know that the rank of the tensor 
defining a four-qutrit state has rank 1269, which 
seems out of reach for any practical computa¬ 
tion. 

• Bell inequalities. 

Bell inequalities related to AME(4,3) are un¬ 
known. 

The small corner of the Hilbert space formed by 
AME states may well be more complex than expected. 
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Appendix A: Four-qubit entangled states 

For completeness we will present in this appendix 
a discussion of maximally entangled states of four 
qubits. To see that no AME(4,2) states exist it is 


sufficient to analyze the purity of reduced density ma¬ 
trices. The purity Trp^, serves as a measure of the 
degree of mixedness of the density matrix p, but also 
as a measure of entanglement of the initially pure state 
reduced to p by a partial trace. 

The same argument works for multipartite systems. 
Let us denote the qubits by A, B, C and D. The purity 
of the bipartition AB of a state \tp) is given by Tr(p^^) 
with pAB = Trc£)|' 0 )(?/;| being the reduced density 
matrix of AB. The theoretical minimum for this quan¬ 
tity is 1/A^; so 1/4 in the case of four qubits. However, 
one cannot attain this minimum for the three biparti¬ 
tions at the same time, as it was proven analytically 
in Ref. [ 6 ]. The best one can have is 1/3 in all bipar¬ 
titions. One state that accomplishes this was given in 
that same paper: 

\HS) = ^ [|0011)- 1 -|1100)-ha;|0101)-f 
v 6 

wllOlO) -ha;2 (10110) -H |1001))] , 

(Al) 

where w = exp(^). This state has also maximum 
entropy of entanglement. Another state that has min¬ 
imum purity in all bipartitions is: 

\HD) = ^(10001)-f 10010)- 1 -10100) 

V 6 

-K|1000)-f V 2 |llll)). (A 2 ) 

This state was found by the authors and by [Hj to 
also carry maximum hyperdeterminant, an extension 
of the concept of determinant to higher dimensions 
|47j . and an interesting entanglement measure that is 
intrinsically multipartite. An equivalent state already 
appeared in as an example of a ’’symmetric maxi¬ 
mally entangled state” of 4 qubits. 

Let us note that AME states are defined through 
their entanglement properties and thus can be trans¬ 
formed into any equivalent state under local unitaries. 
Gour and Wallach [48j found \L) and \M) states 
while searching, respectively, for the states that max¬ 
imize the average Tsallis a-entropy of entanglement 
for a > 2 and for 0 < a < 2 : 
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\L) (((1 + a;)(|0000) + |1111)) + (1 - a;) (|0011) + |1100)) 

+u? (10101) + 10110) + |1001) + |1010))) (A3) 

((^ <!“">+"™»+ 

^ (10101)+ 11010))^, (A4) 


where w = Remarkably, states \L) and \M) can 

be transformed by SLOCC into the states \HD) and 
\HS) respectively. They are called SLOCC equivalent. 

In fact Verstraete et al. gave a classification of all 
pure 4-qubit states in 9 SLOCC inequivalent classes 
|22) . The most important class is called the generic 
class and is presented by Gour and Wallach [iS] in the 
following compact form: 

g = {zo\(j)''~)\(j)''~) + Zi\(l)~)\(j)~) + 

+ Z3\t{j~)\llj~) \ Zo,Zi,Z2,Z3 e C}, (A5) 

where |0^) = (|00) ± |ll)/-\/2 and \^p^) = (|01) ± 
\\Q)/y/h are the Bell states. 

The method to transform a state from the compu¬ 
tational basis into the Bell basis (which means ob¬ 
taining the values of the Zi) is explained in detail 
in Verstraete’s paper. Any 2 states that have the 
same 4 parameters Z 3 ,zi,Z 2 and Z 3 in this basis are 
SLOCC-equivalent. The SLOCC-equivalence between 
the states \L) and \HD) (as well as between the states 
\M) and \HS)) is proven by showing that they have 
precisely the same zi coefficients. 


a vector of length 4 and putting into the 
matrix X^. That is, 


' Xu 

X21 

Ai3 

X23 

X12 

X22 

Xi 4 

A24 

A31 

X42 

X33 

X43 

Lx24 

A42 

X34 

A44 . 

An 

Ai2 

X21 

X22' 

Xi 3 

Xi 4 

X23 

A24 

A31 

X32 

X41 

X42 

_ X33 

X34 

X43 

A44 


reordered 


(Bl) 


Here, colored exchanged entries are set in boldface 
(colors in online version). The other 3 matrices are 
transpositions of U, and U^. Therefore, a matrix 
U of size 4 is multi-unitary if U, and are 
unitary. The same restrictions hold for matrices of 
size (P. 

It is possible to demonstrate that 2-unitary matrices 
of size U = 4 do not exist [5] . The smallest 2-unitary 
matrix exists for order D = 9 - see Eq.(43|. 


Appendix B: Partial transposition, reshuffling 
and 2-unitarity 

For any matrix A of a square order D = dP repre¬ 
sented in a product basis, Xmn ■= (m, n|A|p,, i^), one 

defines [JS] its partial transposition, X^^ = Xmv, and 

reshuffling, = Xm^- To get a better feeling of 

nv 

these particular reorderings of elements we consider 
a matrix X of order four. Let us now switch to the 
standard, two-index notation, and write its elements 
as Xij with = 1,2,3,4. Here, we have ( 2 ) = 6 
different reorderings of 2 indices out of 4. Two of 
these reorderings are particularly interesting: (a) the 
partially transposed matrix is equivalent to the 
matrix with all four blocks of size two transposed, and 
(b) the reshuffled matrix A^ is obtained by taking lex¬ 
icographically each 2x2 block of A, reshaping it into 


Appendix C: Symmetric sudoku, generalized 

AME(4,3) states and 2—unitary permutations 

According to the usual sudoku rules all digits in 
a single row/column of the matrix or block of size 
three are different. Let us distinguish other sets of 
nine elements: location, which contain all digits from 
the same place in each block (e.g. nice centers of the 
blocks - set in red in the matrix below); broken rows 
containing three rows of length three occurring in the 
same position of three blocks (blue example in the 
matrix) and an analogous notion of broken column. 

Note that the standard operation of matrix trans¬ 
pose, A^, exchanges columns of the matrix with its 
rows. The partial transposition A^^ exchanges rows 
with broken columns and columns with broken rows. 
Furthermore, reshuffling operation, A^ interchanges 
blocks with rows and columns with locations - see 
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The following matrix shows an example of a sym¬ 
metric sudoku pattern analyzed in m- all digits in 
each column, each row, each block, each location (red 
in online version), each broken row (in blue), and each 
broken column (boldface) are different. 
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4 
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6 
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(Cl) 


Consider now a matrix Vi of size nine with all en¬ 
tries equal to zero besides nine ele men ts placed in the 
positions of digit ’1’ of the matrix (Cl I equal to unity. 
As all nine digits in each location of Sg are different 
it is clear that this is a legitimate permutation matrix 
of size 9. As the address of each non-zero element can 
be interpreted as a pair of two ternary digits, it repre¬ 
sents a Graeco-Latin square (28) and determines the 
AME(4,3) state |fl) 

Thus permutation matrix Vi is 2-unitary, as its 
partial transpose and reshuffling remain unitary. The 
same property holds also for other permutation ma¬ 
trices Vm with TO = 2,..., 9, obtained by placing nine 
ones in the positions occupied by digits’m’ in pattern 


(Cl). The property of 2-unitarity is preserved if we 
enphase a permutation matrix by multiplying it by 
a diagonal unitary matrix T), and take V!^ — VVm- 
Thus it is fair to say that the symmetric sudoku ma¬ 
trix (Cl) encodes nine families of enphased 2-unitary 
permutation matrices or, families of AME(4,3) states. 


Appendix D: Mini catalogue of multi—unitary 
matrices 

Multi-unitary matrices are defined for orders equal 
to powers of integers. As there are no 2-unitary matri¬ 
ces of size 4 the smallest interesting cases are D = 8,9 
and D = 16. In general, multi-unitary permutation 
matrices are one-to-one connected m to combinato¬ 
rial arrangements called orthogonal arrays of index 


unity [33] ■ Here, we present a resume of all the multi¬ 
unitary matrices presented in this work and some ad¬ 
ditional ones: 


{i) D = 2^ = 8. The orthogonal matrix Og defined 
in (41) is 3-unitary. This matrix is equivalent 
with respect to permutations and enphasing to 
the three-qubit Hadamard matrix iLg = 
(which is 1-unitary). It is not possible to con¬ 
struct a 3-unitary permutation matrix of size 

D = 8 [T2]. 


(ii) D — 8^ = 9. There exist 2-unitary permutation 
matrices related to symmetric sudoku designs - 
see Appendix Also permutation matrices de¬ 
fined at the end of Section [Vl B| and the complex 
Hadamard matrix ( |43| ) are 2-unitary. The latter 
one is equivalent to the tensor product of two 
Fourier matrices Fg 0 Fg, but this product is 
1 -unitary only. 


{in) F = 2^ = 4^ = 16. In this case 

there are no 4-unitary matrices m- There 
exist a 2-unitary permutation matrix given 
by Perm(4, 3,13,10,14, 9, 7, 0,11,12, 2, 5,1, 6,8,15) 
created from OA(16,4,4,2), which was generated 
using the Gendex Module NO A [50] . 


{iv) D = for a prime p. Two-unitary matri¬ 
ces of this size exist, but they are neither com¬ 
plex Hadamard matrices (d^ non-zero entries) 
nor permutations {d non-zero entries). Indeed 
they have p^ non-zero entries which correspond 
to powers of the main root of the unity 
Construction: Let us denote (j, k) to the j-th 
row block of size p and the k-th column block 
of size p oi U (note that U has blocks of size 
p). Each block (j, k) has to be filled with the dis¬ 
placement operator Dj^k of size d. Displacement 
operators are defined in Eq.(46). The explicit U 
for d = 3 is provided in Eq.(451. 


As we can see, for some dimensions there are no 
2-unitary permutation matrices. Indeed, a 2-unitary 
matrix of size exists if and only if a projec¬ 

tive plane of odd order d exist (see Theorem 8.43 in 
|34j 1. In particular, they exist for every prime (and 
odd prime power) d. For example, the existence of 
AME(4,3) states relies on the existence of the projec¬ 
tive plane of order 3. 
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